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TRANSFORMATION OK SURFACE AI.BEDO TO 


SURFACE - ATMOSPHERE ALBEDO AND IRRADIANCE, 

AND THEIR SPECTRAL AND TEMPORAL AVERAGES 

Myron L. Nack 
Robert J. Curran 

ABSTRACT 

The dependence of the albedo at the top of a realistic atmosphere upon the 
surface albedo, solar zenith angle, and cloud optical thickness is examined for 
the cases of clear sky, total cloud cover, and fractional cloud cover. The radi- 
ative transfer calculations of Dave and Braslau (1975) for particular values of 
surface albedo and solar zenith angle, and a single value of cloud optical tl k- 
ness are used as the basis of our parametric albedo model. The question of 
spectral and temporal averages of albedos and reflected irradiances is ad- 
dressed, and unique weighting functions for the spectral and temporal albedo 
averages are developed. 
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TRANSFORMATION OF Sl’RFACE ALBEDO TO 


SURFACE - ATMOSPHERE ALBEDO AND IRRADUNCE, 

AND THEIR SPECTRAL AND TEMPORAL AVERAGES 

1. INTRODUCTION 

Radiation transfer calculations form the coupling between satellite deter- 
mined radiation budget components and the physical processes incorporated in 
climate models. The use of exact radiation transfer calculations in such climate 
models is often frustrated by the volume of tabular and graphical data produced 
in the transfer calculation, and the non-analytic nature of the solution The 
present study is the result of an effort to match the needs of climate modelers to 
the available radiation budget data by means of a parameterization of exact 
radiation transfer calculations. A recent set of calculations of this tvpe was 
performed by Dave and Braslau (1974, 1975), hereafter referred to as DB, for a 
realistic atmosphere with and without a stratus cloud layer. Both sets of calcu- 
lations are here parameterized to relate surface-atmosphere albedo to surface 
albedo in terms of linear equations. 

A linear albedo relationship for the cloud free atmosphere has been developed 
by Vinnikov (1965) and used quite extensively by Budyko (see Budvko (1974) ). 

The linear relationships of Vinnikov are traceable back to 1962; a time when 
measurements of the optical properties of the atmosphere were somewhat in- 
complete. However, that formulation has been found to agree almost exactly 
with the present parameterization for a solar zenith angle 6 O’. In a sense, the 
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present formulation extends the work of Vinnikov to a more general dependence 
on surface albedo, Bolar zenith angle and cloud optical thickness. In sections 
2 and 3 of this paper "cloudy" means total cloud cover and surface albedo is 
assumed to be horizontally homogenous. In section 4 we extend our analysis 
to include fractional cloud covr and fractional surface albedo cover. 

In subsequent work, Curran, Wexler and Nack (1978) t hereafter referred to 
as CWN (1978), apply the derived albedo transformations to a monthly zonal sur- 
face albedo model. This application produces a monthly zonal cloud free and 
cloudy surface-atmosphere system albedo which is used in the determination of 
monthly zonal cloud fractions. The modeled albedos are combined with the mea- 
sured data supplied by the Nimbus-6 Earth Radiation Budget (ERB) experiment 
(see Smith et al., 1977) to form the cloud fractions. These cloud fractions arc 
necessarily determined for a fixed local time due to the sun-synchronous orbit 
of Nimbus-6. 

The final effort in this series Is the determination of the daily and monthly 
averaged, reflected zonal shortwave irradiances. The irradlances are formed 
by assuming the cloud fractions to be constant over al 1 . local times throughout 
the month. Daily integrations of the irradiances are made using the averaging 
techniques developed in this present report. These techniques make use of our 
solar zenith angle parameterization. They could include time dependent cloud 
fraction data, but will use our ERB albedo determined cloud fractions due to a 


lack of time dependent data. 
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2. A LB K DO TRANSFORMATIONS FQH CLEAR AND CLOt’DY 

ATMOSPHERES 

The albedo of the earth-atmosphere system Is defined as the ratio of re- 
flected solar irradtance to that incident from the sun. For the present study 
the system albedo, denoted A s for cloud free and A c for cloudy, is assumed to 
be a function o r the surface albedo a, the solar zenith angle 0 (l and a cloud 
optical thickness r c . The surface albedo a is defined in a manner similar to the 
system albedo, i. e. , the surface albedo a is the ratio of the irradiance reflected 
from a surface to that incident upon it. For the purposes of the present study, 
the surface is assumed to be a Lambert reflector. It should be noted that the 
refectance of natural surfaces may depend on the angular distribution of the 
incident radiation. This latter dependency is assumed to be second order to 
the parameterization being investigated. The cloud optical thickness is defined 
at 0. 555 /im. 


To fit the results of Dave and Braslau (1974, 1975) linear functions are 
chosen as follows; 


A s <a.0 (( ) * 11^(0,, )a ♦ A,(0, 0„>. 


( 1 ) 


A (a.0 n .r c ) 


m t .<0 o )g( T c )a + A t (0 ’°0' V 

A c !a,0 o , A c (0. 0. t )| 


( 2 ) 


Terms dependent on higher powers of a were neglected since the fit resulting 
from the linear functions seemed satisfactory. 

Equation (2) indicates that the explicit r c dependence of A can be replaced bv 
an implicit dependence on \ (0, 0, r c ). The slope factor of Fq. 1 depends 



only on solar zenith angle 0 () . However, the slope factor of I q. 2 has been 
separated into a factor dependent «>n solar zenith angle n\ ( 0 O ), and a factor 
dependent on cloud optical thickness g( r c ). 
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Since the Dave and Braslau results were computed at only one optical thick- 
ness (r c a 3 . 35), the present formulation sets g (3. 35) = 1. 0. To accomodate 
other values of r t or /\ ( 0 , 0 , r ) we require g to decrease as its argument 
increases. A function w'hich satisfies th.s constraint is: 


g(T k ) = g|A to (i. ri| 


\ (U.O, t <5 I 
c 

A (0. U. r l 


(3) 


The purpose of the function g is to dampen the dependence of A c (a, t) 0 , r c ) on a 
as r increases. A further constraint on g is 


V '-W < 10 (4) 

foi all a including the extreme case of a = 1 . 0 , since the realistic atmospheric 
model of Dave and Braslau (1974, 1975) that we parameterize Includes absorption. 

The intercepts A s (0, 0 (| ) and A ( 0 , « (| , 3. 35) of Fqs. 1 and 2 are found using 
Table 4 of Dave and Braslau (1974) for the atmospheric models Cl and Cl -ST, 
respectively. A polynomial fit to these data was made as follows: 

A s «).0 o ) - 0.04k l + ( 1 .OS 7 X 1O" 5 >0* - (2.219 X IO' q )0^ 

♦ ((>.77(i X IO' I3 )0^ 


(5) 


with (0, 0, r. ) * 0. 1940 in this section. The solar zenith angle fl () is here 
measured in degrees. The solar zenith angle independent terms of Eqs. 5 and 
0 are identically equal to the Dav.j and Draslau values for 0 O * 0*. The re- 
maining coefficients of Eqs. 5 and 6 were obtained by performing a least 
squares difference minimization fit of the polyi.omial to the tabular values of 
Dave and Draslau. The solar zenith angle dependent values of A v (0, 0 O ) and 
A c (0, 0„, r t . ) using Eqs. 5 and 6 are shown in Table 1. These intercept values 
are also shown graphically in Fig. 1. 


Table 1 

The Intercepts and Slopes of Linear Albedo Transformations 


°0 

A, 


m s 

n\ 

0 

0. 04 83 

0. 1940 

0. 7213 

0. 5079 

5 

0. 0480 

0. 1952 

0. 7213 

0. 5037 

10 

0. 0494 

0. 1989 

0. 7213 

0. 4977 

15 

0. 0500 

0. 2050 

0. 7212 

0. 4900 

20 

0. 0523 

0. 21.% 

0. 7209 

0. 4804 

25 

0. 0544 

0. 2247 

0. 7203 

0. 4091 

30 

0. 0508 

0. 2382 

0. 7192 

0. 4500 

35 

0. 0595 

0. 2541 

0. 7171 

0. 4411 
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Table 1 (Continued) 


°0 

_~-T 


1 

i 

in, 

n\ 

40 

0. 06 28 

0. 2725 


0. 7137 

0. 4245 

45 

0. 066 8 

0 2933 


0. 7083 

0. 4060 

50 

0. 0722 

0. 3167 


0. 7000 

0 3858 

55 

0. 0796 

0. 3424 


0. 6877 

0. 36 38 

60 

0. 0903 

0. 3706 


0. 6700 

0. 3400 

65 

0. 1057 

0. 4013 


0. 6451 

0. 3144 

70 

0. 1280 

0 4 344 


0.6108 

0. 2870 

75 

0. 1598 

0. 4 700 


0. 5644 

0. 2579 

80 

0. 204G 

0 5080 


0. 5025 

0. 2270 

85 

0. 2666 

0. 5485 


0. 4212 

0. 1943 

90 

0. 3509 

0. 5914 


0. 3157 

0. 1598 


The slope functions m s and n\ of I qs. 1 and 2 are determined by using the 
above intercept values and making a best fit to the surface albedo dependent results 
of Dave and Braslau. The best fit is for the angles 0 (J = 0“ and 80° of Table 7 
of Dave and Braslau (1 97*4 ), and 0 O = GO" in Fig. 7 of Dave and Braslau (1975). 

The values obtained by this fitting procedure are shown in Table 1 for these 
three angles. 

The functional relationships which were used to fit the 0 O dependence of 

the slope functions for these three angles are giv<*n bv: 

m s (0 o ) = 0.7213 - (2.180 X IO' 9 >0f, - <4.‘>4I X 10“' 3 . (7) 

m c (0 o ) = 0.5079 - (h.59(i X IO' 4 »0 (| - (3.5t>5 X IO' 5 )0j- . 


(8) 








ALBEDO 
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Figure 1. Solar Zenith Angle Dependence of Our Intercept Albedos 
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The three coefficients of in* (0„) ami m*. ( 0 O ) are determined by fitting their 
valueu at 0*. 80\ and 80* exactly. Graphs of the slope functions of Kqs. 7 and 
8 are shown In Fig. 2, anil their values are shown In Table 1. 

The 0 ,, dependence of both the slopes and Intercepts of Kqs. 1 and 2 may 
now be substituted Into those expressions to form the top of the atmosphere 
albedo for both the cloudy and cloud-free atmosphere. The values of A s (a. 0 O ), 
the cloud-free albedo, are given In Table 2 as a function of both arguments. 

These values are also given In graphical form In Fig, 3. Similar values for the 
albedo at the top of the cloud obscured atmosphere are shown In Table 3 and 
displayed graphically In Figure 4. 

If we compare the values of A, (a, 0„ ) In Table 2 with the values of A c (a, 

0„ . r , ) of Table 3 we see that In general A c (a, 0„ , r e ) > A, (a, 0„ ). This 
agrees with our initial intuitive response that the Insertion of the cloud layer 
should increase the albedo at the top of the atmosphere over the corresponding 
lcar sky case. However, a detailed comparison shows that the values above 
the discriminator line in Table 3 satisfies this inequality, whereas, for values 
below the discriminator line this inequality becomes reversed and A s (a, 0„) > 

A c (a, 0 (1 , r c ). This is caused by the absorption In the atmospheric model. The 
optical thickness of the cloud is small enough to allow an amount of radiation 
through the cloud (for smaller solar zenith angles) to be repeatedly reflected 
back and forth between the highly reflective surface and the cloud bottom. These 
multiple reflections in the absorbing atmosphere decrease the amount of 
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radiation that can escape out of the top of atmosphere relative to the clear sky 
case. When 0„ Is large enough, for each of these values of large surface reflec- 
tivity, eventually A c >A,. This Is because the effective optical thickness of the 
cloud gets large enough to make negligible the proportion of radiation penetrating 
the cloud und feeling the effect of these multiple reflections with absorption. 

The separation of values of (a, 0 O ) in Tables 2 and 3 into sets where 
A c > A, and A, > A c leaves a line of values (as shown in Table 3) which satisfy 
A, (a, 0 O ) = A c (a, 0„ t r t ). For these values the ability to distinguish a snow 
covered surface with or without our cloud layer over it disappears if only albedo 
measurements are used. 

It is pci -ji In Figures 3 and 4 that the top of the atmosphere albedo is always 
leas than unity. In agreement with the fact that the radiation transfer calculations 
upon which these pnrameterizations were made, contain contributions due to 
molecular absorption. The increase in the intercepts with increasing Holar 
zenith angle, or dA, (0, 0,, )/d0 o > 0 and OA, (0, 0„ , r t . )/30 o > 0, is shown In 
Figs. 1,3, and 4, and explained by the increase in scattered radiance due to the 
greater effective optical depth ol' an atmosphere in which scattering dominates 
absorption. As 0 () increases less radiation reaches the surface. This addi- 
tional radiation is not lost as a result of the zero surface albedo, and conse- 
quently, is available to the atmosphere for either absorption or scattering. For 
the same reason the slopes decrease as the solar zenith angle increases as is 
shown in Fig. 2, or 3m, (0 O )/ <*0 O < 0 and Dm t (0 (| )/d0 o < 0. This is consistent 
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with the sensitivites of the albedos at the top of the atmosphere to the surface al- 
bedo, or 

3A t (.,0 o )/3 a ■ m f (0 o ) 

3A c (a,0 0 ,r c )/3a • i» c (0 o )g<r c ), 

decreasing us less radiation reaches the surface (due to increasing solar zenith 
angle). 

In general for a fixed surface albedo the albedos at the top of the atmosphere 
increase as the solar zenith angle increases, or dA^a, 0 () )/dd (| > 0 and 
3 A c (a, 0 O , r c )/ 30 o > 0, and this is shown in Figs. 3 and 4. The one exception 
to this relationship is shown in Fig. 3 and occurs for a > 0.8 and 0 Q > 40° when 
A, (a, 0 O ) decreases as 0„ increases. This effect begins to occur for 0 O > 50° 
and a ■ 0.65 as is shown in ’’'able 2. This is because given the two competing 
effects of A 4 (0, 0J increasing and decreasing as 0 () increases, the latter 

effect dominates for large enough a and 0 (l to cause 

< 0 

for this small range of a and 0 () . In this case, the large surface albedo (a > 0.8) 
is a major cause of the large albedo at the top of the atmosphere A, (a, 0 o ) so 
that as 0 n increases beyond 40° less radiation reaches the surrace to take ad- 
vantage of its effect. This results in the albedo at the top of the atmosphere 
decreasing as the effect of the surface albedo decreases with increasing o t) . 


P M W 

300 L d0 0 ' 


3 A (0, 
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In Fig. 3 it was observed by Jacobowitz (1978) that for a surface albedo of 
approximately 0.75 the dependence of A t (a, 0 O ) oo 0 (( is minimal. This feature 
could be used in calibrating aircraft or satellite instruments over areas such as 
White Sands, New Mexico. For example, if we had a three satellite system 
measuring albedo in the morning, noon, and afternoon as proposed by Woerner 
(1978) for the Karth Radiation Budget Satellite System (ERBSS) this independence 
on 0 (| would allow good data comparison for the cloud free data. 

Equation 1 has the property that it may be Inverted to give the surface 
albedo In terms of the albedo at the top of the atmosphere: 

a «V,,'«o' “ lAmeaJV " K<V (9) 

Thus, given a cloud free icene and the measured albedo A, ne> , (0„ ), one can 
solve for the spatially averaged surface albedo in the field of view of the meas- 
urement by using Eqs. 5 and 7 to compute the intercept and slope functions, re- 
spectively. In principle, Eq. 2 could also be inverted to supply surface albedo 
for cloud obscured fields of view. In practice, except for moderate to high 
resolution radiometers, it is difficult to find fields of view which are uniform in 
their atmospheric properties. As a result of variability’ in the physical char- 
acteristics of atmospheric aerosols, such as optical thickness, scattering phase 
function and single scatter albedo, Eq. 9 does not exactly model all atmospheric 
variability accurately enough to permit formation of detailed determinations of 
surface albedo. However, with a relaxation of the accuracy requirements for 
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Figure 2. Solar Zenith Angle Dependence of Slope Functions 



Albedo at the Top of the Cloud Free Atmosphere as a Function of Surface Albedo, 

a, and Solar Zenith Angle, 0 
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Table 3 (Continued) 
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Figure 4. Albedo at the Top of a Cloudy Atmosphere as a Function of 
Surface Albedo and Solar Zenith Angle 
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the determination of surface albedo, global maps of surface albedo can be 
formed using Kq. 9 and satellite determined, high resolution aiix dos. 

3 TRANS FORMATION FOR G ENERALIZED CLOUDY ATMOSPHERE 

As noted in Section 2 both the slope and intercept functions of Kq. 2 should 
depend on the cloud layer optical depth by being explicitly dependent on r c or 
implicitly dependent on r c through a dependence on A t (0, 0, r c ). We in 
general use the latter dependency. 

The intercept dependency of \ (0, 0 {) , r c ) on A c (0, 0, ) is shown in 

Eq. 6. Uslnj; this equation we compute values for A (0, 0 () , t c ) for input 
values of V (0, 0, r c ) = 0.1940, 0.2091, 0.30, 0.40, 0.45, 0.50, 0.55, nnd 
0.00, and these results are shown in Table 4 and Figure 5, The values of 
A t (0, %, r c ) computed for A t (0, 0, r c ) = 0. 2091 are in close agreement 
with the Dl-ST results of Table 4 from DR (1974), even though the increase 
in optical depth in this atmospheric model is due to the aerosol distribution 
and not the cloud layer. 

The dependency of the slope function on t c is contained in the g (r c ) function 
of Eq. 3. 

We now use this implicit dependency on r c of A (a, fl Q , u ) through 
A t . (0, 0, t c ) to compute A k (a, 0 O , t c ) for different values of a and 0 () . 

These results are shown in Table 5 and Figure 6 for A t (0, 0, r.) = 0. 50. 

We again see that the albedo at the top of the atmosphere has a very weak 
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Table 4 

Cloud Albedos, \ (0, 0 () , r. ) 
\ (0, 0, rj 


% 

0. 19 

0. 21 

0. 30 

0. 40 

0. 45 

0. 50 

0. 55 

0.60 

0 

0. 19 

0. 21 

0. 30 

0. 40 
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0. 50 

0. 55 

0.60 
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0. 19 
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Figure 5. Cloud Albedos as a Function of Solar Zenith Angle for the Five 
Different Intercept Values A (0, 0, r c ) of Table 4. 




Albedo, A (a, <? 0 , r. ), at the Top of The Cloudy Atmosphere as a Function of Surface Albedo, a and 
Solar Zenith Angle, 0 V with the Cloud Optical Thickness, r , , Contained in the Input Parameter 

A. (0, 0, t ) = 0. 5. 
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dependence on surface albedo for large solar zenith angles; e. g. A,., (a, 90, r t ) 
increases only 0.79 f when a increases from 0.0 to 1.0 We also see that in gen- 
eral a comparison between Tables 5 and 2 phows that A t . > A, , except for values 
below the discriminator line in Table 5 when A, (a, 0„ ) >A t (a, 0„ , r e ). The 
reasons for this effect are discussed in Section 2. We see that the range of 
values of a and 0 o when A, > A c are now decreased with an increase in the cloud 
optical thickness, when A c (0, 0, r t ) increases from 0.194 in Table 3 to 0.50 in 
Table 5. 

4. SPATIAL RESOLUTION AND FRACTIONAL SURFACE 
AND CLOUI) COVER 

A general measurement of albedo at the top of the atmosphere occurs 
during conditions when a fraction of the field of view (FOV) of the instrument 
is covered by clouds (f ) and a fraction of the surface in the FOV is covered by 

C 

an albedo a k , where the fractions (f k ) of the k = 1, N types of surfaces satisfy 

N 

1 " £ f k' (1°) 

k=l 

We will assume that the surface albedo in the FOV can be expressed by 

N 

k=l 


a 


(ID 




2(3 


The application of our cloud free (A 4 ) and total cloud cover (A ) albedos to the 
analysis of albedo data (A) with fractional cloud cover will now be presented for 
the cases of high and low spatial resolution albedo data. We will assume we 
have this data at time t, with a solar zenith angle 0 Q (t), over an average surface 
albedo a obtained from Kq. 11. The cases we discuss are high and low resolu- 
tion data, fractional cloud cover with random and known distributions over the 

FOV, and data for these values of a, 0 U , and t averaged over a time period T. 

If we have high spatial resolution data so that the surface in the FOV can be 
characterized by a single value of surface albedo and the size of the FOV allows a 
yes or no cloud decision then the albedo A for a given £icture element (pixel) can 
be used to classify that pixel using the following tests; 

Test Clas s of Pixe l 

A < A 4 (a, 0 (l ) cloud free, low aerosols 

A, (a, 0 o ) < A < A t . (a, 0„, r c ) cloud free, heavy aerosols, 

or small clouds 

A c (a, 0 O , r c ) < A cloud 

where A ( (0, 0, r ) =0.20, for example, could be used to compute A (a, n , r ) 
for our lower cloud threshold in this test. The cloud pixels could be further di- 
vided into low, medium, or high albedo (or optical thickness) clouds if additional 

tests on A using selected ranges of (0.20, 0.35), (0.35, 0.5C), and >0.50 
for A t . (0, 0, r. ) were used to compute test ranges of A (a, 0 (| , t c ). 
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If we have low spatial resolution data where clouds of characteristic 
optical thickness r t . or albedo A (0, 0, T Mn general fill a fraction f of the 
FOV then the explicit albedo dependence on our variables can be expressed by 

A - A(a,fl 0 .r c ,f c ) (12) 

A simplified way of representing the dependence of A on (a, 0 ()I r ) implicitly 
is given by 


A * A(A s .A c .f c ) * A|A 1 (a,f 0 ).A c (a,f 0 .r e ) l f c |. 


(13) 


If we assume the form of Eq. 13 is valid then a polynomial in f of the form 


a I 

n c 


(14) 


“n “ «n<\ A c> = Vl 


(15) 


is one way of approximating A. If the value of A (0, 0, r t . ) selected to charac- 
terize our typical cloud causes A (a, 0 () , r. ) to be an upper limit to the measured 
values of A then three constraints which could be used to determine « n are 


A = \ if 

f t = o, 

(16) 

A = A c if 

f. = 1. 

C 

(17) 

A, < A < A c if 

0 < f. < 1. 

C 

(18) 


The Monte Carlo results of Davis, Cox, and McKee (1978) for finite and in- 
finite cloud layers support the decrease in A of Eq. 18 for f c < 1. In Fig. 7 
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we schematically show the two endpoints of A determined by Eqs. IB and 17, 
ami a linear and two quadratic examples of Eq. 14, For the linear case of 
Eq. 14 the constrainlJ of Eqs. 16 and 17 determine the solution 


% 3 \* 


(19) 


a, « A c - 


( 20 ) 


A = A + ( A - A If = (I - I \)A ♦ I A 

% i* % t C ' t C l 


( 21 ) 



Figure 7. Schematic Relationships Between Measured 
Albedo and Cloud Fraction 


This equation is identical to Fq. 2. 90 of Budvko (1974), and similar to Fqs. 12, 

17, 18, and 19 of Nack and C.reen (1974) who characterized fluxes reaching the 
ground in this manner. For f v = 1/2 in Fq. 21 we have A = 1/2 (A^ + A c ). 

If some technique could determine the value of A for f ( = 1/2, or A (A , A v . , 1 /2), 
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then this value could be used to determine the following solution to the quadratic 
case of Eq. 14; 

®o * V (22) 

a, - 4A(A t . -i\ - A c . , (23) 

a, - 2(A^ ♦ \) - 4A(A % . A^.fc). , (24) 

A = A s ♦ |4A(A s . A c .Vi) - - AJ f, 

♦ 2|A, ♦ A e - 2A(A,,A C ,H)| fj. . (25) 

When A (A,, A c , 1/2 (A, + A c ) this quadratic relationship degenerates 
to the linear relationship of Eq. 21. 

We note that the constraints of Eqs. 16 and 17 are always valid, but that for 
a small range of (a, 0„) the relationship A s > A c can hold as is discussed in 
Sections 2 and 3. When this occurs we must interchange A s with A c in Eqs. 18 to 
25 and Fig. 7. Since the work of CWN(78) uses values of A c close to those of 
Table 5 and values of a < 0. 75 the problem of A, > A c will not arise, j .id will 
not be discussed further. 

A technique which could be used to compute values of A for different values 
of f , such as A (A , A , 1/2), is the Monte Carlo computer simulation of 
radiative transfer through an atmosphere with fractional cloud cover. These 
values of A for different values cf f could then be used to determine higher 
order values of a n . Due to the lengthy computer time requirements of such 
a program it might initially be run for just one characteristic wavelength and 



30 

for a single typical set of values of a, 0 O , and r t using a single random cloud 
distribution over the FOV. 


If auxiliary information was available which gave the spatial distribution 
of the clouds over the FOV then the appropriate surface albedos should be used 
in Kq. 13. For example, if the surface in the FOV is half ice and half water 
with albedos a, = 0. 7 and a tt = 0. 1 for the ice and water, then in general the 
albedos for the cases (1) clouds only over ice, (2) clouds only over water, or 
(3) clouds randomly distributed over ice and water for f c - 1/2 arc all unequal; 


A, |AJa w . 0 „ (. A,. (a,. 0 O . T t . I. ';| * A, | AJa,. 0„>. A (a* . 0„. - c ). *-| 
i A 3 |A % (a. 0 (| ». A t (a. 0 O . T t ). 'i| 


(2G) 


where a = 1/2 (a, + a^ ) = 0. 4 using Kq. 11. 

If at a given time both low and high spatial resolution albedo data were 
available then the fraction of high resolution pixels classified as clouds of all 
the pixels in the low resolution FOV could be used as the truth value of f c . 

This truth value could be compared with the determination of f c from the low 
resolution data found by solving, for example, Kq, 21 for f c , or 

f. = (A - A,)/(A C - (27) 

If we have low resolution albedo data averaged over a time period T 
for the same values of (a, ) then the time averaged albedo data A could be 



31 


expressed tn terms of the time averaged fractional cloud cover f using Eq. 11 
in general, or for example Eq. 21 as follows; 

A|A,.A c .7 c <V) ■ I' - \ + WV (28) 

The temporal averaging of high resolution albedo data Involves direct averaging 
of the results of the cloud or no cloud decision for a given pixel over the period T. 

5. ALBEDO DEFINITION. REFLECTED IRRADIANCE, AND 

SPECT RAL AND T EMPORAL AVI- RACKS 

The albedo of the earth-atmosphere system A is the ratio of the reflected 
irradiance W R over the incident solar irradiance W 0 at the top of the atmos- 
sphere, or A = W K AV 0 . In general, the irradiance quantities are of greater 
significance to studies of the energetics of the atmosphere thar. albedo. This 
definition will be applied to the spectral plane albedo A (X, 0 (( ), and will 
determine the weighting furetions used to define the spectrally averaged albedo 
A (« 0 ), and the temporally averaged albedo A, where ,\ is a particular wave- 
length in the shortwave part of the spectrum used by Dave and Braslau (1974, 
1975) and bounded by X, = 0. 285 ym and X, = 2. 5 *im. The solar zenith angle 
0 () depends on the latitude, the local time t elapsed since noon, and the day of 
the year. Temporal averaging in general would apply U the 0 o (t) dependence 
over a sunlit part of the day during an interval T = t, - t, bounded by (tj , t, ), 
and when these bounds are sunrise and sunset a daily average is obtained. The 
non-weighted spectral < A (0„)> and temporal <A> albedo averages will also 
be given, but these do not satisfy our definition of albedo. 
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The spectrally and temporally integrated reflected irradiances are defined 


as 





W R (X,* 0 )dX, 


(29) 




K 




w K uydi. 


The incident solar radiance 


(30) 


l(X.p.0:g n ,0 o ) ■ H 0 <\> &l/i - - 0„) (31) 

traveling in direction (a, 0) with the sun at direction (/i () , $ 0 ) is lelatcd to the 
incident solar irradiance by 


W 0 (X.fl„) * d.u y* H 0 <X) &(/i - a o )A<0 - 0„) 
Hq H(,(X), 


(32) 


where 6 denotes the Dirac delta function. A direction is specified by (a, <>), 
where u = cos 0 and (0, $) are (zenith, azimuth) angles. The solar energy 
dK 0 incident on an area do over the wavelength and time inter valo d> and dt 
is given by 

dt 0 

"" ,M ■ ' W ° (M > (33) 


The spectrally and temporally integrated incident solar irradiances are given by 


r x 2 

W OS ,0 O> J W n (X.0 o )dX = h 0 H k 


( 34 ) 
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M 


os 



2 

H (| <X)dX - tX 2 


X| t ii )(S 


(35) 


W 


OS 


(«I - 


V 


•'ll 


os'^o ^ 1 


^0 **os 


(3fi) 


" «2 ' i » d ' (37) 

where H os is the incident shortwave energy over intervals dt and do, H ()S is 
the mean value of H os , W os is the incident shortwave irradiance over the 
Interval (t, - t, ), and /J 0 is the mean value of ^ over this time interval. 

To derive our weighting functions we start with our definition applied to the 
spectral albedo 


AtX.tf 


V x V V x -V 
w „< X V ’ «o‘V x> 


(38) 


The weighting function which satisfies our definition for the spectrally averaged 
albedo 


W rx 2 rw„(Mo >1 w r(Mo 

\ I I / I ■ I 

0 w «( 0 o ) A, L wj w 0 o > 

I /* X 2 p , 0 <X, "l 

■ (X, - x.i/ 


dX 


A(X, (t () t dX 


(39) 
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is given by H,, (X)/H 0 . , and this is consistent with Eq. 2 of Dave and Breslau 
(1975). The weighting function which satisfies our definition for the temporally 
averaged albedo 


A 




w 

w 


dl 



A(fl 0 )dl 


(40) 


is given by n 0 /i* n . 

.. s 

The relationships between irradiances and albedcte are now 

W R (\.0 w i = p,, I l M < X » A(X,0 o ), 

/ 

If non-wcighted spectrally and temporally averaged albedos defined by 


(41) 

(42) 

(43) 
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art* used then they have no simple direct relationship with their corresponding 
lr rad lances as In Fqs. 42 and 43. 

We can now combine Fqs. 1 through 8 with Fq. 42 to obtain the following 
transformations of surface albedo to reflected Irradiances at the top of dear 
\V S and cloudy W c atmospheres; 

w s < a *V " V. (46) 

W " "o l, 0sV‘®o r c>* ( 47 ) 

If we apply Fq. 21 to the fractional cloud cover case we have 

W(a. # 0 .T C , f e ) - g () H s A(A J .A c .f c ) 


■ (I - f c > W s ta. ® 0 > ♦ 




(46) 


* W(W.. W . I I. 

set 

where in general we expect 0 () , r t , and f to depend on the time of day t at a 
given location. 

6. DAILY AND MONTHLY AVERAGED ALBEDOS AND REKLKCTED 
IHRAPIANCKS 

One of the basic problems in acquiring climatological radiation data from 
satellite platforms is the temporal sampling problem. A polar orbiting satel- 
lite will sample the reflected irradiance which is used to compute albedos at a 



3(j 


given time of day for many locations so that global spatial sampling is achieved. 
A geostationary satellite could simple « given region of the globe with high 
temporal resolution, so that a costly system of well spaced geostationary satel- 
lites could give fine temporal and spatial sampling resolution. We will now show 
how our 0 o (t) parameteri/ations of albedos can be used to overcome the tem- 
poral sampling problem of the less costly nolar orbiting satellite. 

We wish to compute time averages of albedo and reflected irradiance at the 
top of the atmosphere at a latitude and longitude for an area whose surface 
albedo, a m) , remains constant during the J 1 *' nay of the m 11 ' month. The cloud 
free and total cloud cover cases will be presented fir°t where r c will be in- 
dependent of time. The time dependence of the solar zenith angle 0,,(t) can be 
computed from a well known equation of astronomy Isee p. 50 of Smart (19G5)| 


0„U) ■ 

cor 1 *i„(ti 

(49) 

u,,(t) = sin6 m) sin^ 

♦ cns6 m | cos0 cos lit t * 

(50) 

lilt) 

* 2 tr t/T 

(51) 


where 0 is the latitude, 6 m) is the solar declination, T is the earth’s rotation 
period (24 hours), t is the time elapsed since noon, and H is the hour angle. 
The times of sunrise and sunset, ±T niJ , are defined by p,,(t) = 0 so that Eq. 

50 gives us 

cos(2»T ni j/T) = - tan6 m jtan0. 


(52) 



and the minimum solar zenith angle 0 m) Is determined by t » 0 so that 

cos0 m j ■ sin6 mJ wn0 ♦ cosA m jVOt# • costO - A m j). 

°m) m <t> - * n ij • 

If we let (t , , t : ) - (-T nii , T niJ ) In Kq. 37 we can obtain the dally averaged 
values of the cloud free A s (a m| , t mj ) and total cloud cover A^ (a mJ , t n)) , t ) 
albedos by using Kq. 40 as follows; 

/ “ t in j 

/i„(t)dl. (65) 

> mi 

ll l Mu 1 1 *P|H| I dt « (60) 

>n,) 

/ •'mi 

l i *nij , ®o(^ , ^c I !(*«(* * ^*in) I dt , (->7) 

l,„j 

where t m) indicates the value of a fun' ion averaged over the j' 11 day of the in 1,1 
month, and t m the value of the function averaged over the j - 1 . . . . , N m days 

of the m lh month which has N m days. Equation 43 can now be used to relate 

% ‘ 

these daily averaged albedos with daily averaged irradtances resulting In 


( 53 ) 

(64) 
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where W os (t mj ) 1* the daily averaged incident solar irradiance of Eq. 30, and 
the dally variation of !l os is noted due to the elliptical orbit of the eprth about the 
sun causing the radius distance to vary. 


in general the monthly averaged reflected and incident solar irradiances 
are given by 


W R (, m> * rr“£ 

' i i 


(61) 


'm 


m j.i 


(62) 


The weighting function w'hich satisfies our definition for the monthly averaged 
albedo 


A(t m ) = 





N m 




^mj Hos u mj * \ 


(63) 


is therefore given by 


w n:j ^ in j ^OS^mj^^mj ^*os^mj^ • 


(04) 


If we apply Eqs. 61 to 64 to the cloud free and total cloud cover cases w’e obtain 
the following monthly averaged irradiances and albedos: 


w s(*m 1 

^rnj ^kis^mj ^ ) 

(65) 

W c (t m ) 

^mj ^os (*mj ^ * *iii ^ 

(66) 
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As^ni* “ v. w mj A<, ^ a mj'*mj * • 

"in j-l 


(G7) 


A^'m* “ vi ^ w mj ^c^mj 1 *mj ,r c ^ • 
'« j-l 


( 68 ) 


We now address the problem of temporal averaging when fractional cloud 
cover occurs throughout the day. The dally averaged lrradlancc is obtained by 
combining Fqs. 21, 30, and 48 resulting in 

1 




ox 

••■mi 


f III - f c (OI W s (a mj ,0 o ) + f c (t)W c (a mj ,0 o ,r c (t))\ dt(69) 


m J *'- T >ni ^ 


where in general f c (t) and r c (t) have some diurnal dependence on t. if we use 
values of f c and r c which are constant over the day, but have mean values over 
the month, f c (t m ) and r c (t m ), which may vary from month to month then 


I * “ )) W s (a m j,t m j) + f c (t m )W C l a mj'*inj’ T C (*m 

^mj^os^mj^ ^ ^ As(a tn j,l m j ) + l c (t m )A C |a ni j,t m j.r c (t, n )|^ 

and the monthly averaged ir radiance is 


(70) 


W(t m ) = M - fc<t m )iw s (t m ) + f c Ct m )W c (t m ) 

^ ^mj^os (*mj ^ ^ )l As (t m ) ■*" ( c (tm )A C I tpi ,r c^m » 

where the daily and monthly averaged albedos are 


(71) 
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A|amj.t m j,T c (tn 1 ).tc(tni)l ~ 1 1 - f c (t m )| A s (a m j,t m j) 

♦ f c (t m )A c |a m j,t m j,r c (t m )| 
A(l m ) - 1 1 - f c (t m )| A s (t m ) + l c (i m )A c |t m ,r c (l m )| . 


(72) 


(73) 


in CWN (1978) we use surface albedo climatology data, our parameterizations 
of A s (a, 0„ ) and A t (a, 0„, r c ), and the albedo model of Eq. 21 for total albedo, 
A(a, 0„, r c , f c ), for an atmosphere with fractional cloud cover, f t , to compute 
monthly and zonally averaged values of f c . Our values of A are obtained from 
the Nimbus 6 Earth Radiation Budget (ERB) instrument which is discussed by 
Smith et al. , (1977). Consequently, our values of f c are only representative of 
the Nimbus G sampling times. If the diurnal behavior of cloud cover was Unown 
then it could be fit to our values of f t . at these times to give f t . (t), and Eq. G9 
could be used to compute the daily averaged irradiance. Since this diurnal 
behavior is not known, in a subsequent paper we will use the values of f c (t m ) 
determined in CWN (1978) in Eqs. 70 to 73 to determine daily and monthly 
averaged values of irradiance and albedo. 

Since Eq. 27 was used to determine f c at the ERB sampling time we know 
that the value of albedo and irradiance are consistent with the ERB values at 
that time of day. Our use of our parameterizations of A s (a, 0 „ ) and A c (a, 0 o , 
t c ) at other times of the day in Eqs. 56, 57, and 70 through 73 allows us to 
approximately solve the temporal sampling limitations of Nimbus G and compute 
proper time averages of albedo and irradiance. 
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Due to the simplicity of our 0 () parameterizattons of A s (a, 0 O ) and A c . (a, 

0 O , r c ) we will only need to compute #J m j and the weighted mean values of 0"(t) 
defined by 

<Cj = TT“ / 0*>l*><‘)/*mjldt (74) 

- r mj J- T mJ 

for n » i, 2, 4, and (5 in order to obtain the daily averaged values of Eqs. 56 
and 57. If we view the time dependent albedos to be specifically dependent on 
0" (t) as a polynomial Instead of generally dependent on 0 n (t) then the daily 
averaged values of A s |a mj , 0"(t)| and A c |a mj , 0"(t), t c | are obtained by sub- 
stituting 0[J,j for 0"(t) in these functions, or 

A c ( a m j.t m jl ~ \ia m j,0 m j,r t ) . (76) 


7. CONCLUSIONS 

We have developed a simple parameterization of the extensive radiative 
transfer calculations of Dave and Braslau (1974, 1975) which should prove use- 
ful to many scientists in the areas of earth radiation budget, climate modeling, 
and climatological data analysis, and we ourselves will use these results in our 
future work. We have noted that the single cloud model of DB needs generali- 
zation and have attempted a reasonable approach to this problem. 
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The temporal and spatial sampling problem of satellite experiments designed 
to measure time averaged albedos and irradiances at the top of the atmosphere 
is discussed. We present how the use of our solar zenith angle parameter! zatlons 
can allow proper time averages to be performed given satellite albedo measure- 
ments sampled at a few discrete times. By applying the definition of albedo 
consistently we derive the weighting functions needed to perform proper spectral 
and temporal averages of albedo. 
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